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Abstract 

We explore the influence of in-medium nucleon-nucleon cross section, symmetry potential and 
impact parameter on isospin sensitive observables in intermediate-energy heavy-ion collisions with 
the ImQMD05 code, a modified version of Quantum Molecular Dynamics model. At incident 
velocities above the Fermi velocity, we find that the density dependence of symmetry potential plays 
a more important role on the double neutron to proton ratio DR(n/p) and the isospin transport 
ratio Ri than the in-medium nucleon-nucleon cross sections, provided that the latter are constrained 
to a fixed total NN collision rate. We also explore both DR(n/p) and Ri as a function of the impact 
parameter. Since the copious production of intermediate mass fragments is a distinguishing feature 
of intermediate-energy heavy-ion collisions, we examine the isospin transport ratios constructed 
from different groups of fragments. We find that the values of the isospin transport ratios for 
projectile rapidity fragments with Z > 20 are greater than those constructed from the entire 
projectile rapidity source. We believe experimental investigations of this phenomenon can be 
performed. These may provide significant tests of fragmentation time scales predicted by ImQMD 
calculations. 
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I. INTRODUCTION 



The nuclear symmetry energy plays an important role in the properties of nuclei and 
neutron stars To a good approximantion, it can be written as 

E sym = S(p)5 2 . (1) 

where 5 = (p n — p p )/(p n + P P ), is the isospin asymmetry; p n , p p , are the neutron, proton 
densities, and S(p) describes the density dependence of the symmetry energy. Theoretical 
predictions for S(p) from microscopic nucleon-nucleon interactions show large uncertainties, 

n q 

especially in the region of suprasaturation density [6|, |7[ • Constraining the density depen- 



dence of the symmetry energy has become one of the main 



goals in nuc 



27H35, 



41 



ear physics and 



42( | . Heavy Ion 



has stimulated many theoretical and experimental studies [7H19I. 
Collisions (HIC) with asymmetric nuclei provide a unique opportunity for laboratory studies 
of the density dependence of the symmetry energy because a large range of densities can be 
momentarily achieved during HICs. The experimental strategy is to measure the collisions 
of neutron-rich and neutron-poor system and compare them in order to create observable 
that are primarily sensitive to the density dependence of the symmetry energy. Comparisons 
of these observables to theoretical predictions provide the sought-after constraints. 

Several semi-classical transport theoretical models have been developed to to simulate the 
nucleus-nucleus collisions and have been used to obtain constraints on the symmetry energy. 
In Section II. A, we discuss and compare the approaches using the Boltzmann-Uehling- 
Uhlenbeck (BUU) equation and the Molecular Dynamics Model (QMD). In this work, we 
choose to simulate nuclear collisions with the code ImQMD05 developed at the Chinese 



Institute of Atomic Energy (CIAE) 42l444j . Details of this code are described in section II. 
B. In Section II. C, we show the dynamics of the simulated collisions of U2,12A Sn -|- 112 > 124 Sn 
at E/A = 50MeV. The influences of the density dependence of symmetry energy, in- medium 
nucleon-nucleon cross sections, and impact parameters on isospin sensitive observables for 
the simulated reactions are investigated in Section III. We focus our studies mainly on 
two experimental observables. We examine the ratio between neutron and proton yields 
in Section III. A, and the isospin transport ratio in Section III.B where we also discuss the 
influence of cluster emission on isospin transport ratio. Since the main goal of this paper 
is to investigate the influences of various parameters used as input for the transport model 



there is no attempt to obtain the best fit to the data as was done in Ref . 33] . We summarize 
our findings in section IV. 



II. TRANSPORT DESCRIPTION OF HEAVY ION COLLISIONS 

One frequently utilized transport models to describe the heavy ion collisions is the 
Boltzmann-Uehling-Uhlenbeck (BUU) equation, which provides an approximate Wigner 
transform of the one-body density matrix as its solution [36(. Another frequently utilized 
approaches, known as the Molecular Dynamics Model (QMD) represent the individual nu- 
cleons as Gaussian "wave-packet" with mean values that move in according the Ehrenfest 
theorem; i.e. Hamilton's equations [45]. Even though we do not use the BUU equations in 
our work, it is useful to discuss the similarities and differences between the two approaches. 



A. Remarks on QMD and BUU 

At the code level, both BUU and QMD models propagate particles classically under 
the influence of a mean field potential, which is calculated self-consistently the positions 
and momenta of the particles, and allow scattering by nucleon-nucleon collisions due to 
the residual interaction. The Pauli principle in both approaches is enforced by application 
of Pauli blocking factors. These similarities in implementation have lead to similarities in 
predictions for many collision observables |[37j. 

There are also significant differences in these approaches. In the BUU equations, each 
nucleon is represented by 200-1000 test particels that generate the mean field and suffer 
the collisions. In QMD, there is one test particle per nucleon. A-body correlations and 
cluster formation are not native to the original BUU approach; which is supposed to provide 
the Wigner transform of the one body density matrix. On the other hand, many-body 
correlations and fluctuations can arise from the A-body dynamics of QMD approach. Such 
A-body correlations are suppressed in BUU approach, but correlations can arise in both 



approaches from the amplification of mean field instabilities in spinodal region [38]. Collision 
algorithms in the QMD approach modify the momenta of individual nucleons, while in BUU 
approach, only the momenta of test particles are modified. Depending on the details of 
the in-medium cross sections that are implemented, the blocking of collisions can also be 



more restrictive for QMD than for BUU, leading to fewer collisions and therefore a greater 
transparency. Since there are typically more than 100 test particles per nucleon, collision 
induced fluctuations are smaller in BUU than in QMD possibly suppressing the fragment 
formation rates. 

Fragments can be formed in QMD approaches due to the A-body correlations and these 
correlations are mapped onto the asymptotic final fragments by a spanning tree algorithm. 
Serval different methods have been developed to allow BUU codes to calculate cluster pro- 
duction. In the Stochastic Mean Field (SMF) approach, for example, the time evolution of 
the one-body phase-space distribution / is governed by the nuclear mean-field, two-body 



scattering, and a fluctuating (stochastic) term, which causes the fragmentation 2114241] . Such 
fluctuations are sometimes achieved by reducing the number of test particles in a BUU sim- 
ulation to about 50 per nucleon. The pBUU code developed in Ref. [26] provides a sophis- 
ticated theoretical description of clusters up to A=3, but lacks heavier particles, including 
alphas. Also, the coalescence approximation or the spanning tree method has predicted 
light cluster yields from the sampled test particle distributions. Of these three methods, the 
production of clusters has a feedback on the collision dynamics in the first two methods, 
while the coalescence or the spanning tree methods do not. 

Without modifications to produce clusters, direct comparisons of BUU calculations to 
the observables for Intermediate Mass Fragments (IMF's), such as the rapidity dependence 
of isospin transport ratios extracted from the yields of mirror nuclear fragments, are rarely 
performed. Such fragment observables have been modeled within Stochastic Mean Fields 
(SMF) models, which add fluctuations to the BUU and within the QMD approach. 

There are many versions of BUU and QMD codes used by different groups. Even though 
there is substantial overlap between the results, the extracted constraints on the symmetry 



energy based on isospin diffusion data obtained from QMD approach (ImQMD05 (4 lN44| ) . 



tend to favor a somewhat softer density dependence of symmetry energy term in the equation 



of state than those current 
SMF 
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20|, 



y provided using BUU approaches, such as the IBUU04 19 
22(, and pBUU(26| codes. 
In some case, the BUU approaches predict much smaller values for the double n/p ratios, 
DR(n/p), than the experimental results. It is important to determine to what extent these 
differences arise from (1) differences in the underlying transport models like those discussed 
in the preceding paragraphs, or from (2) differences in the mean fields or in medium cross 



sections used in these calculations or (3) in their numerical implementation in the codes. 

Differences in the results between the various QMD and BUU approaches can arise from 
the transport parameters used in the codes. The two main ingredients in transport mod- 
els, the nucleonic mean fields and nucleon-nucleon binary scattering cross sections, are not 
treated consistently with the same microscopic nucleon-nucleon interactions used for both 
mean field and collisions. Thus, there can be differences in the approximations made to 
adapt these quantities to the calculations and associated differences in the predictions made 
with them. While the effects of the symmetry energy around saturation density and the 
isospin dependence of in-medium nucleon-nucleon cross sections have been studied in several 
BUU codes 19], it has not been similarly studied in QMD models. This paper provides a con- 
strained exploration on the sensitivity to the isospin dependence of in-medium corrections 



by expanding upon the ImQMD05 calculations published in Ref. 



33|. 



B. Brief description of ImQMD05 



Early implementations of the QMD model did not explore the density dependence of nu- 
clear symmetry potential jssj. At CIAE, we developed and successfully applied a new QMD 
code, labeled ImQMD, to study heavy ion reactions close to the Coulomb barrier. Such as 
the excitation function for fusion cross sections |40j. In this code, we incorporated a symme- 
try energy potential that depended linearly on density. Modifications incorporated in the 
later ImQMD05 version of this code include mean field potentials calculated using a Skyrme 
energy density functional with options for different forms of the density dependence of the 



symmetry potential 
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44j . With these modifications, ImQMD05 has successfully described 



the multiplicity of reaction products, collective flows and stopping powers in intermediate 
energy Heavy Ion Collisions (HICs). More recently, it has been used to constrain the density 
dependence of symmetry energy at sub-saturation density using three experimental observ- 
ables: (1) double n/p ratios, (2) isospin diffusion constructed with isoscaling parameters, 
and (3) with the ratio of mirror nuclei as a function of rapidity for 112 > 124 Sn -f 112 > 124 Sn at 
E/A=50MeV[33j. 

In this section, we will describe in more detail the theoretical formalisms used in the 
ImQMD05 code 4jJ-|44j. Within this code, the nucleonic mean fields acting on nucleon 



wavepackets are derived from a potential energy density functional where the potential 
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energy U includes the full Skyrme potential energy with the spin-orbit term omitted: 

U = U P + U md + U coul (2) 

Here, U cou i is the Coulomb energy, while the nuclear contributions can be represented in 
local form with 

U Pi md = J u P!md d 3 r (3) 

and, 

a P . ft P^ + 9sur /r7 n2 

u p =77— + i -i fj + 77— Vp 
2 p 7? + 1 Po 2 Po 



Po 



: (V(p„ -p p )) 



2 



n8/3 



2 Po Pq 7 

where the asymmetry 5 = (p n — p p )/(p n + p p ), and p n , p p are the neutron and proton 
densities. In the present work, a symmetry potential energy density of the form (p^) 7 ^ 2 P 
is used in transport model comparisons. The energy density associated with the mean-field 
momentum dependence is represented by 

1 ^ 1 



I I r 

Umd= — T7J~6 / d3 Pl d3 P2fN 1 (Pl)fN 2 (P2) 



2 Po Artlv 2 167r 

1.57[ln(l + 5 x l(T 4 (Ap) 2 )] 2 . (5) 

Here, Jn are nucleon Wigner functions, Ap = |pi — P2I , and the energy is given in MeV 
and the momenta in MeV/c. The resulting interaction between wavepackets in Eq. (J5]) 



is described in Ref. 45j. All of the calculations use the potential given above, and the 
coefficients used in formula (Hj) are directly related to the standard Skyrme interaction 
parameters, with a = — 356MeV, (3 = 303MeV, r] — 7/6, g sur = 19A7MeV fm 2 , g sur ,iso = 
-11.35MeVfm 2 , C s = 35.19MeV, and g pT = 0. Eq. (JSJ modifies the EoS so that it 
no longer matches the EoS obtained with the Skyrme interaction parameters at T=0MeV. 
Appendix discuss the influence on isospin observabels when the isoscalar part of the EOS 
obtained with Skyrme interaction is changed by the introduction of momentum dependence 
to Eq. (jSJ). Within statistical error, we found that it has no significant effect on the isospin 
observables studied in this paper. 

The isospin-dependent in-medium nucleon-nucleon scattering cross sections in the colli- 
sion term are assumed to be the form: cr* m/np = (1 - £,{E beam )p/ po)a^ np , where £,{E beam ) = 



0.2 for Eb eam = 50AMeV, and the isospin dependent o~^f np is taken from Ref. 46]. The 
in-medium isospin dependent nucleon-nucleon (NN) differential cross section in free space is 



46| , and isospin dependent Pauli blocking effects are the same as in 
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adopted from Ref. 

n 

44J . Asymptotic clusters are constructed by means of the minimum tree spanning method 
widely used in QMD calculations, in which particles with relative momentum smaller than 
Po and relative distance smaller than i?o coalesce into one cluster. In the present work, the 
values of i? — 3.5/m and P = 250MeV/c are employed. 

From the adopted interaction in the ImQMD05, we construct the density dependence of 
symmetry energy for cold nuclear matter as follows: 

where m is the nucleon mass and the symmetry coefficient C s = 35.19MeV . Using this 
particular parameterization, the symmetry energy at subsaturation densities increases with 
decreasing ji, while the opposite is true for supranormal densities. In general, the EoS is 
labeled as stiff-asy for 7, > 1, and as soft-asy for 7^ < 1. In the present work, calculations 
are performed using 7^ = 0.5 and 2.0 as representative cases of a soft asy-EOS and a stiff 
asy-EoS, respectively. For illustration, the density dependencies of symmetry energy for 
7i = 0.5(solid line) and 7, = 2.0(dashed line) are plotted in Fig. 1. For this work, the region 
of interest is mainly at subsaturation densities, p < 0.16/m -3 . 



C. Reaction Dynamics at E/A = 50 MeV 

Before discussing the influence of the symmetry potential and the in-medium NN cross 
sections on various HIC observables, we first show some predictions for dynamical processes 
and fragment production mechanisms in Sn + Sn collisions at E/A = 50MeV. Figs. 2 (a) and 
(b) show the time evolution of the density contour plots for 124 Sn + 124 Sn at E/A = 50MeV 
for b = (top panels) and b = 6fm (bottom panels) from one typical event. Here, 7, = 0.5 
was used in the calculations, but the density contour plots are practically the same for 
all 7i > 0.2. (For smaller 7$, many more nucleons are emitted during simulations due to 
the instability of the initial nuclei.) The projectile and target touch around 50 fm/c and 
nucleons start to transfer between the two nuclei. The compressed region reaches the highest 
density at approximately 100 fm/c for both impact parameters. However, higher densities 

8 




FIG. 1. (Color online) The symmetry energy for cold nuclear matter plotted as a function of density 
for 7j = 0.5 and 2.0 (solid line) for symmetry potentials used in ImQMD05 simulations (Eq. [6]). 
For reference, we show the symmetry energies used in IBUU04, for x=0 and -1 (dot dashed lines) 
from Ref [29] . 

are reached in the overlap region and more compressional energy is stored for calculations 
with the soft-asy EoS potentials than with the stiff-asy EoS potentials. 

In central collisions (6 = Ofm) at an incident energy of E/A = 50MeV, the system 
reaches a maximum density at about lOOfm/c and subsequently expands to low density 
where multifragmentation occurs. At this point, the reaction system with the soft-asy EoS 
disintegrates into more light fragments than the system with the stiff-asy EoS. In peripheral 
collisions (6 = 6/m), the projectile and target residues separate after 150/m/c has elapsed, 
whereupon the low-density neck connecting them ruptures into fragments. The two excited 
residues then continue along their paths without further mutual interactions. Thus the 
nucleon diffusion process in peripheral heavy ion collision is predicted to terminate at about 
150/m/c. 

Fig. 3 shows the charge (left panels) and mass (right panels) distributions of the reaction 
products for 112 Sn + 112 Sn at an incident energy of 50MeV per nucleon for 6 = 2,6, 8fm. 
The solid symbols are the results for 7, = 2.0 and open circles are the results for 7, = 0.5. 
The charge and mass distributions depend on the symmetry potential parameter 7^. The 
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FIG. 2. (Color online) Time evolution of the nucleon density in the reaction plane for 124 Sn + 124 Sn 
collisions at E/A=50MeV at b=0 (top panels) and 6 (bottom panels) fm. 



distributions obtained with 7, = 0.5 are narrower than those obtained with 7^ = 2.0. 

For completeness, the multiplicities of fragments with charge Z as a function of their 
scaled rapidity y/yl'J£ m for 112 Sn + 112 Sn and b = 2,6,8/m are shown in Fig. 4. The 
maximum charge of fragments at mid-rapidity decreases with increasing impact parameter b. 
In general, two ridges of heavier fragments are observed, distributed near the initial projectile 
and target rapidities. For reference, the initial projectile and target rapidities are marked by 
the vertical dotted lines at y/yt'^ m = 1 and —1, respectively. The heaviest fragments have 
lost about 35% of their initial velocity for central collisions, and about 25%(10%) of their 
initial velocity for b = 6fm (b = 8fm). The velocity loss of the heaviest fragments depends 
on the decelerating effects from the effective N-N interactions and nucleon-nucleon collision 
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471 ] . Intermediate mass 



frequency. These features have been observed in experiments 
fragments with Z between 3 and 20 are produced over a wide range of rapidities for every 
impact parameter. 



III. ISOSPIN OBSERVABLES 

In this section, we examine the uncertainties in constraining the density dependence of 
the symmetry energy from transport model simulations. To do this, we study the effect 
of varying the isospin dependence of in-medium NN cross-sections, the impact parameter, 
and the symmetry potential on isospin sensitive observables. In particular, we focus on the 
calculated results of the neutron to proton yield ratio as well as the isospin transport ratio 
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FIG. 3. (Color online) The charge (left panels) and mass (right panels) distributions for 
112 Sn + 112 Sn at b=2 (top panels), 6 (middle panels), and 8 (bottom panels) fm. The solid symbols 
are for ji = 2.0 and open symbols are for ji = 0.5. 




beam 



FIG. 4. (Color online) Multiplicity of fragments with charge Z as a function of their scaled rapidity 
V/Vblam for 112 Sn+ 112 Sn at E/A=50MeV for b=2 fm (left panel), 6 fm (middle panel), and 8 fm 
(right panel). 



(Ri). These two observables have been used extensively to 
dependence of symmetry energy at subsaturation density 



provide 
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constraints on the density 



23 



28|. 



A. Nucleon yield ratios 



The symmetry potential, which is opposite in sign for protons and neutrons, directly 



affects the emission o 
quite intensely 



28 



29 



nuc 
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eons. This fact was recognized early on and has been studied 
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52j. Experimentally, due to their low interaction with mat- 
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ter, neutrons are more difficult to measure. Furthermore, different types of detectors are 
usually employed to measure neutrons and protons. To minimize the effects of detector 
efficiencies, the yield ratio from an asymmetric system A is compared to the yields from a 
more symmetric system B by constructing the double ratio DR(n/p) = R n / P (A) / R n / p (B) = 
(Y n (A)/Y p (A))/(Y n (B)/Y p (B)), where Y n (A) and Y P (A) are the neutron and proton yields 
obtained in system A, and Y n (B) and Y P (B) are the neutron and proton yields obtained 
in system B. In Fig. 5, the DR(n/p) data from two collision systems A = 124 Sn + 124 Sn 
and B = 112 Sn + 112 Sn for central collisions are plotted as solid stars in the left panel. The 
shaded regions in Fig. 5 are the ImQMD05 results for two different symmetry energy cases 
7i = 0.5 (upper shaded region) and 2.0 (lower shaded region), where the transverse emitted 



4l|. 



nucleons have been selected with an angular cut of 70° < 9 c . m , < 110° as in Ref. 

The DR(n/p) values for the 7, = 0.5 case are greater than that for 7^ = 2.0. In inter- 
mediate energy HICs, the emitted nucleons mainly come from the overlap region during the 
expansion phase. This region is below saturation density, where the symmetry energy is 
larger for smaller 73. Larger symmetry energy in this region results in enhanced neutron 
emissions from the neutron-rich system and thus larger values of the double ratio. This 



behavior has been observed in other transport model predictions [281. 1291. 152j]. However, some 
of these other calculations using BUU approach produce double ratios that are much smaller 
than the data for all symmetry energies, while the ImQMD05 calculations for 73 = 0.5 are in 
the correct range. It remains a puzzling discrepancy between transport model predictions, 
and detailed comparisons between the different codes and different transport models should 
be done in the future to understand this difference better. Such study would allow more 
definite constraints on the symmetry energy from transport models. 

The impact parameter dependence of DR(n/p) ratios is plotted in the right panel of 
Fig. 5. The energy cut, 40 < E cm ^ < 80MeV, and angular gates are chosen to minimize 
sequential decays and cluster effects on the calculated values. The cuts also provide more 
robust coalescence-invariant quantities. The impact parameter dependence is rather weak. 
(Noting that there is a large zero-offset on the y-axis.) One reason for the insensitivity with 
respect to impact parameter in these calculations is that transverse emitted nucleons mainly 
come from the overlap region, where the high energy values of DR(n/p) are related to the 
information of symmetry energy of nuclear matter, regardless of the size of that region. 

Similar to the symmetry potential, the isospin dependence of in-medium nucleon-nucleon 
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FIG. 5. (Color online) (Left panel) Double n/p yield ratios, DR(n/p), from transverse emitted 
nucleons as a function of kinetic energy. The solid stars are data from Ref. [2^]. (Right panel) 
DR(n/p) from transverse emitted high kinetic energy nucleons (-E7 c .m. > AQMeV) as a function of 
impact parameter. The upper shaded region is for 7$ = 0.5 and the lower shaded region is for 
7i = 2.0. 

(NN) cross sections should play an important role on the isospin sensitive observables for 
heavy ion collisions in the transport model. The in-medium NN cross sections depend on 
the type of colliding particle, the relative energy of the colliding pair and the local density of 
medium. The calculations are performed at b=2fm where the sensitivity to the in-medium 
cross section are relatively larger because the density is higher on the average and there are 
more nucleon-nucleon collisions for central nucleus-nucleus collisions. We investigate three 
options for the in- medium NN cross sections :1) we used the free nucleon-nucleon cross 
sections, with no reduction cr^ n ^ np = cr(J^L p ', 2) We used an experimentally motivated a-^ N 
reduction, cr* n / np = (1 — 0.2 p/p )a^ e np ; 3) We used isospin independent reduced cr* n/ , np case, 
<n = °; P = < P = (1 - 0.2p/poK, where a' = (2N np a^ + (N nn + N pp )a f J n % p )/N NN . In 
case 3, N np , N nn / pp , N NN are the number of possible combinations of np, nn, pp, and NN 
colliding pairs in the reaction system. By keeping the total number of NN collisions to be 
almost the same when we compare calculations with case 2 to those with case 3, we can 
study the effect of the isospin dependence of the in-medium cross sections without varying 
the overall nucleon-nucleon collision rate. Here, we are explicitly interested in the isospin 
dependence. 

The calculated results of the DR(n/p) ratios are plotted in Fig. 6 for the energy and 

13 



angular gates described above, which emphasize the emission from the overlap region formed 
during the early stages of the HICs. By comparing case 2 (open squares) to case 3 (solid 
circles), we can test whether it makes any difference to the double ratios, whether the overall 
collision rate is dominated by n-p collisions as one would expect if the in-medium reduction 
is the same for all cross sections, or whether n-n, p-p and n-p collisions are governed by the 
same in-medium cross sections. The effect of reducing the cross sections below their free 
values can be investigated by comparing the DR(n/p) results obtained with case 1 (inverted 
triangles) and case 2 (open circles). The changes in the cross section have no effect for 
ji—0.5. For 7j = 2.0, there is a slight difference between case 1 and case 2, indicating 
some influence of the total number of collisions when the symmetry potential is smaller at 
subsaturation density. However, the effect of the isospin dependence of the in-medium NN 
cross sections on the DR(n/p) is smaller than the statistical errors for both 7$ we studied. 
Overall, the DR(n/p) exhibits little to no dependence on the in-medium NN cross section. 
The main reason is that NN collisions are suppressed by Pauli blocking at E/A = 50MeV. 
The different &nn/np/ a nli/np or a n P / a nn values used in the transport model do not lead to 
significant differences in the DR(n/p) ratios. Thus, the isospin observable DR(n/p) from 
Sn+Sn at E/A = 50MeV is a robust observable to constrain the density dependence of 
symmetry energy. 



B. Isospin transport ratio R, t 

When the projectile and target nuclei come into contact, there can be exchange of nucleons 
between them. If the neutron to proton ratios of the projectile and target differ greatly, the 
net nucleon flux can cause a diffusion of the asymmetry 5 reducing the difference between 
the asymmetries of the two nuclei. This isospin diffusion process, which depends on the 
magnitude of the symmetry energy, affects the isospin asymmetry of the projectile anc. 
target residues in peripheral HICs. The isospin transport ratio Ri has been introduced [l4] 
to quantify the isospin diffusion effects, 

Ri = a \ b \ (7) 

where X is an isospin observable and the subscripts a and b represent the neutron rich and 
neutron-poor nuclei. In this work, we use a and b to denote the projectile (first index) and 
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FIG. 6. (Color online) Effects of different in-medium NN cross sections on the DR(n/p) from 
transverse emitted high energy nucleons (-Ec.m. > 40MeU). The inverted triangles are for the 
free space NN cross section (case 1 in text). The open squares are for the isospin-dependent in- 
medium NN cross section with a phenomenological formula (case 2). The solid circles are for the 
isospin- independent in-medium NN cross section (case 3). 



target (second index) combination, where a = 124 Sn, and b = 112 Sn. We obtain the value of 
Ri by comparing three reaction systems, a + a, b + b and a + b (or b + a). Construction of 
the transport ratio minimize the influence of other effects besides isospin diffusion effects on 
the fragment yields, such as preequilibrium emission and secondary decay, by rescaling the 
observable X for the asymmetric a+b system by its values for the neutron-rich and neutron- 
deficient symmetric systems, which do not experience isospin diffusion. Based on Eq. ([7j), 
one expects Ri = ±1 in the absence of isospin diffusion and Ri ~ if isospin equilibrium is 
achieved. Eq. (J7|) also dictates that two different observables, X, will give the same results 
if they are linearly related. In one experiment, X was taken as the isoscaling parameter, a, 
obtained from the yie.d of the light partic.es near the projeeti.e rapidityfl. In transport 
models |14j . |15| based on the Boltzmann Uehling Uhlenbeck (BUU) approach, the isospin 
asymmetry 5 of the projectile residues (emitting source) has been used to compute Ri(S). 



The constraints on the density dependence of symmetr y en ergy were obtained by comparing 
BUU predictions to the isospin diffusion data 19l . |23| . This is possible only if the 



isoscaling parameters a are linearly related to the isospin asymmetry, 5, o 



residues. This has been shown both experimentally and theoretically 
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collisions. 

In the BUU analysis, the projectile residues are defined as the sum of regions with local 
density greater than 0.05po and center of mass velocity greater than half of the beam velocity 
[14J]. The definition of a residue based on this criterion is not applicable to the QMD model 
as a clustering process is introduced and these residues disintegrate, producing fragments. In 
our ImQMD05 calculations at b = 6fm and E/A = 50MeV, we find that isospin diffusion for 
124 Sn + 112 Sn at b = Qfm stops around 150/m/c. The emission of nucleons and fragments 
at later stages is not part of the isospin diffusion process, but may change the isospin 
contents of the final fragments and thus needs to be included in the analysis in order to 
accurately quantify the diffusion process. In the following, we calculate the collisions of 
four reaction system, 124 Sn + 124 Sn, 124 Sn + 112 Sn, 112 Sn + 124 Sn, 112 Sn + 112 Sn at incident 
energy of 50MeV per nucleon for variety of parameters. We use two different isospin tracers 
to calculate the isospin transport ratios of Eq. ([7]). Our studies lead us to propose another 
isospin transport ratio for constraining the symmetry energy in future experiments. 

We first analyze the amount of isospin diffusion by constructing a tracer from the isospin 
asymmetry of all emitted nucleons (N) and fragments (frag), including the heavy residue if it 
exists, with velocity cut v^^ raa > 0.5v£ e ™ m (nearly identical results are obtained with higher 
velocity cut v^^ rag > 0.7v^ m ). This represents the full projectile-like emitting source, 
and should be comparable to what has been measured in experiments. Fig. 7 shows the 
results of isospin transport ratios Ri(X = SNjrag) (upright triangles) as a function of the 
impact parameter for a soft symmetry case (7^ = 0.5, open symbols) and a stiff symmetry 
case (7j = 2.0, closed symbols). Ri obtained with soft-symmetry case is smaller than those 
obtained with stiff-symmetry potential case. This is consistent with the expectation that 
higher symmetry energy at subnormal density leads to larger isospin diffusion effects (smaller 
Ri values). 

Ri depends weakly on impact parameter over a range extending from central (b = 3/m) 
to mid peripheral^ = 8fm) collisions. Interestingly, the isospin equilibrium and global 
thermal equilibrium are not reached even for central collisions. It is different than the re- 



sults obtained in 23]. Our results show, that neither the effective interaction is sufficiently 
strong nor the collisions are sufficiently frequent (most of them are Pauli suppressed) to 
mix the projectile and target nucleons completely. These two effects prevent the combined 
system from attaining isospin equilibrium even in central collisions. With impact parameter 
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FIG. 7. (Color online) Isospin transport ratios as a function of impact parameter with two tracers 
for a soft symmetry case (7$ = 0.5, open symbols) and a stiff symmetry case (7$ = 2.0, closed 
symbols). Upright triangle symbols are for the tracer defined by the isospin asymmetry of all 
fragments and unbound nucleons with velocity cut [v^^ ra9 > 0.5t;£ e ™ m ), X = 5j^j rag . Circles are 
for the tracer defined by the heaviest fragment with Z max > 20 in projectile region, X = 5z max >20- 

increasing for b > 5fm, the overlap region and thus the number of nucleons transferred 
from projectile and target decreases, causing the Ri values to increase. The impact param- 
eter dependence of Ri predicted by ImQMD05 has been compared to experimental data at 
35AMeV where the measured trends are more consistent with the softer symmetry energy 
(7i = 0.5) than with the stiffer symmetry energy (7$ = 2.0) {5]]]. 




In peripheral collisions, most often, a large residue remains. If it decouples from the full 
emitting source before it equilibrates, it may experience a different amount of diffusion than 
the full emitting source examined by X = d^jrag- To examine this, we constructed a tracer 
using the isospin asymmetry of the heaviest fragments with charge Z max > 20 in the pro- 
jectile region. This tracer is mainly relevant to peripheral collisions as the central collisions 
are dominated by multifragmentation and very few large projectile fragments survive. The 
dependence of Ri(X = 5z max >2o) for impact parameter b > 5fm is shown as open and closed 
circles in Fig. 7. The isospin transport ratios constructed from the different isospin tracers 
have different values especially in the case of 7^ = 0.5. Stronger isospin equilibration (smaller 
Ri values) is observed in the isospin transport ratios Ri(X = S^jrag) constructed from nu- 
cleons and fragments than Ri(X = Sz max >2o) constructed from the heaviest fragments with 
Zmax > 20. Since isospin diffusion mainly occurs through the low-density neck region, and 
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the system breaks up before isospin equilibrium, the asymmetry of the projectile and target 
residues do not achieve equilibrium and, larger Ri(X = 5z max >2o) values result. In contrast, 
there is more mixing of nucleons from the target and projectile in the neck region due to 
the isospin diffusion. Consequently, rupture of the neutron-rich neck is predicted to result 
in the production of neutron-rich fragments at mid rapidity. Experiments [53( have been 
proposed in both Rare Isotope Beam Factory at RIKEN and the Coupled Cyclotron Facility 
at NSCL to test these predictions by comparing the isospin transport ratios obtained from 
residues and from intermediate mass fragments at different rapidity values. 

Since fragments are formed at all rapidities, we can examine the rapidity dependence of 
Ri to obtain more information about the reaction dynamics. Fig. 8 shows Ri as a function 
of the scaled rapidity y/ybeam- The symbols in the leftmost panel are experimental data 
obtained in Ref. [l8] for three centrality gates. This transport ratio was generated using the 
isospin tracer X = ln(Y( 7 Li) / Y ( 7 Be) ) where Y( 7 Li)/Y( 7 Be) is the yield ratio of the mirror 



nuclei, 7 Li and 7 Be (18J. As expected the values of Ri obtained from peripheral collisions 
(solid stars) are larger than those obtained in central collisions (open stars). For comparison, 
the ImQMD05 calculations of Ri(X = S^jrag) are plotted as lines in the middle and right 
panels for a range of impact parameters. The middle panel contains the results from the 
soft symmetry potential (7$ = 0.5) while the right panel shows the results from the stiff 
symmetry potential (7$ = 2.0). The impact parameter trends and magnitude of the data 
are more similar to the results of the calculations from soft symmetry potentials (7, = 0.5), 



consistent with previous analysis 33(. However, the experimental trend of Ri gated on the 
most central collisions (open stars) is not reproduced by the calculations. The experimental 
data indicate more equilibration for central collisions near mid rapidity while the transport 
model indicates more transparency. The equilibration in the E/A = 50MeV data may be the 
result of the impact parameter determination from charged particle multiplicity wherein the 
most central collisions are assumed to be the ones with highest charge particle multiplicity. 
For the most central events, a gate on the highest multiplicity, may select events in which 
more nucleon-nucleon collisions occur rather than a strict selection on the most central 
impact parameters. We plan to study the impact parameter smearing effect in detail in a 
future paper. 

Next, we discuss the influence of the isospin dependence of the in- medium nucleon-nucleon 
cross section on isospin diffusion, using the three cases described in Section III. A. Previous 
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FIG. 8. (Color online) (Left panel) Experimental Ri as a function of rapidity for three centrality 
gates [16]. (Middle panel) The calculated results of Ri(X = 8 n, frag) as a function of rapidity for 
6 = 2,4, 6, 8/ to for 7$ = 0.5 and (Right panel) 7, = 2.0. 



studies have shown that larger in- medium NN cross sections enhance isospin diffusion 48], 
and so the value of Ri in case 1 (free cross sections, inverted triangles) should be smaller 
than the Ri values for case 2 (reduced cross sections, open squares). But as Fig. 9 shows, 
there is no noticeable effect at the energy we studied. Two effects lead to this result. First, 
the dynamical process of heavy ion collisions is governed more by the mean field than the 
nucleon-nucleon collisions at 50MeV per nucleon. Second, the nucleon-nucleon collision 
frequency is small for peripheral HICs. The direct comparison between the results from 
BUU and ImQMD is not appropriate as the two calculations use very different density 
dependence of symmetry energy and in-medium nucleon-nucleon cross section. In fact when 
similar range of symmetry stiffness is used in both IBUU04 (x=-l and x=0) and ImQMD05 
(7i=0.5 and 2.0) as shown in Figure 1, the magnitude of the Ri values are similar. Any 
remaining differences could be caused by the differences on the exact form of in-medium 
NN cross sections and Pauli blocking adopted in the simulations. Indeed, both transport 
models predict that isospin diffusion depends more strongly on the symmetry potential and 
less strongly on the isospin dependence of the in-medium NN cross section at this beam 
energy. 
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FIG. 9. (Color online) The effects of different in-medium NN cross sections on Ri(X = 5^j rag ). 
The inverted triangles are for the free space NN cross section v^n/rip = G^njup ( case 1 i n text)- The 
open squares are for the isospin-dependent in-medium NN cross section with phenomenological 
formula o"* n/ / np = (1 — 0.2p/ po)a^ e np (case 2). The solid circles are for the isospin-independent 
in-medium NN cross section o m ^ in — ®pp — ® 

* (case 3). 



IV. SUMMARY 



In summary, we have investigated the influences of the density dependence of the sym- 
metry energy, of the in-medium nucleon-nucleon cross section and of the impact parameters 
on several different isospin observables predicted by ImQMD05. The study shows that these 
isospin dependent observables are more strongly influenced by the mean field than by the 
NN collisions. In particular, the double n/p yield ratios (DR(n/p)) and the isospin trans- 
port ratios (Ri) are more sensitive to the density dependence of the symmetry energy than 
to the isospin dependence of the in-medium NN cross sections. This conclusion is similar 
to conclusions reached using BUU approaches in the range of symmetry energies studied 
here. The fact that two different approaches used to simulate isospin diffusion draw the 
same conclusions indicates real progress on constraining the symmetry energy with heavy 
ion collisions. 

We also examine the impact parameter dependences of the DR(n/p) ratio and of the 
isospin transport ratio Ri. Our results illustrate that there is only a weak impact parameter 
dependence of DR(n/p) from central to peripheral collisions. Likewise, there is only a weak 
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impact parameter dependence of Ri from central to mid-peripheral collisions. From mid- 
peripheral to peripheral collisions, however, the predicted Ri values increase rapidly with 
impact parameter reflecting the decreasing number of nucleons transferred between projectile 
and target in increasingly peripheral collisions. 

Cluster formation is important for intermediate energy heavy ion collisions. The 
ImQMD05 approach produces clusters at all rapidities, which allows us to study reaction 
dynamics by means of Ri as a function of rapidity. We find that the ImQMD05 approach de- 
scribes the data better with soft symmetry energies term than with stiff symmetry energies. 
We also tested different tracers by constructing corresponding isospin transport ratios for 
them using different symmetry energies. For weakly density dependent symmetry energies 
(small 7j) with large symmetry energies at sub-saturation densities, the values of Ri for the 
residue tracer X = Sz max >20 are larger than those extracted from the entire emitting source, 
i.e., X = 5]yj rag . Calculations using the latter tracer has been previously compared to the 
measured isospin diffusion data in heavy ion collisions. The difference between these two 
tracers can be examined experimentally as a new probe of the symmetry energy. 
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Appendix: I 

The explicit momentum dependent interaction (MDI) term introduced in Eq. (j5]) con- 
tributes to the EoS at T=0MeV as follows: 

Emd/A= Emdl P = 

^[18( M V + l)logV/ 3 e+l) 
36e d 

-6(2u 2 e 3 - 3w 4/3 e 2 + 6u 2/3 e + 11) log(w 2/3 e + 1)) 

+ (4u 2 e 3 - 15u 4/3 e 2 + QQu 2/3 e)} (A.l) 
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Here u = pj po, e = 21.57 and the energy is in MeV. This form of the momentum dependent 
interaction changes the mean field potential and therefore the EoS. Eq. (JSJ) should then 
have a correction to recover the EoS obtained with the Skyrme interaction parameters at 
T=0MeV, which is not included in our calculations in this paper. To see the influence of 
this correction, we can add a counter term as follows. 

1 ^ 1 

1.57[ln(l + 5 x l(T 4 (Ap) 2 )] 2 



1 1 f 

Umd = — TeT~6 / d3 Pl d3 P2fN 1 (Pl)fN 2 (P2) 



"° N lt N 2 

1.57[ln(l + 5 x l(T 4 (Ap) 2 )] 2 (A.2) 



Here, /at (pi) are the nucleon Wigner functions at a given temperature, and the f^ipi) 
represent the nucleon Wigner density at T=0MeV. The second term is just the first term 
evaluated at zero temperature, so with this correction to the momentum dependent inter- 
action we recover the EoS of cold nuclear matter (T=0MeV) obtained with the Skyrme 
interaction parameters. 

When we repeat the simulations of the Sn + Sn collisions using the momentum dependent 
interaction of Eq. ( 1A.2I) . albeit with less statistics, the results of the mass/charge distribu- 
tions, neutron/proton ratios and isospin transport ratios are the same as the results obtained 
using Eq. (jHJ) without the correction term, within statistical uncertainties. Therefore, the 
two different momentum dependent interactions in Eq. Q and Eq. flA.2j) do not lead to 
significant differences in the observables. 
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